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THEORETICAL  PARACHUTE  INVESTIGATIONS 
Progress  Report  No  21 

INTRODUCTION 

1.0  This  Is  the  twenty -first  quarterly  report  covering 

the  time  from  1  March  1962  to  31  May  1962  on  the  study  program 
on  basic  Information  of  Aerodynamic  Deceleration. 

1.1  As  In  preceding  reporting  periods,  work  during  this 

reporting  period  has  been  pursued  In  accordance  with  the 
technical  program,  and  is  described  in  the  following  sections 
of  this  report. 
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Project  No  1 


2.0 


2.1 


Investigation  of  Wake  Effects  on  the  Behavior  of 
Parachutes  and  Other  Retardation  Devices  Behind 
Large  Bodies  at  Subsonic  and  Supersonic  Speeds 


Width  of  Wake  Behind  Basic  Bodies  of  Revolution  in 
Subsonic  Plow 


2.1.1  Introduction 

Progress  Report  No  20  presented  initial  results  of  an 
attempt  to  determine  experimentally  the  wake  width  close  behind 
ten  bodies  of  revolution  in  subsonic  flow.  It  was  postulated 
that  a  functional  dependence  of  wake  width  on  distance  down¬ 
stream  could  be  expressed  in  the  form 


where  b  «  radial  width  of  wake 

R  *  radius  of  primary  body  ■=  D/2 
n,k  *=  constants 
X  =  distance  behind  body. 

The  constants  in  this  empirical  relationship  were  determined 
for  each  body  from  experimental  data  presented  initially  in 
Progress  Report  No  17. 

In  this  report,  the  experimental  velocity  distribu¬ 
tions  for  the  ten  bodies  are  presented  more  accurately,  and 
the  details  of  the  derivation  of  the  empirical  constants  are 
presented.  In  addition,  similar  information  is  presented  for 
an  ogive  cylinder  with  no  afterbody,  based  on  experimental 
data  obtained  from  Progress  Report  No  17*  and  also  for  a 
circular  disk  based  on  experimental  data  from  Ref  1. 
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2.1.2 


Analysis 


Tables  1-1  through  1-12  present  experimental  velocity 
distributions  In  the  wakes  of  the  twelve  bodies  of  revolution 
shown  in  Pig  1-1  for  X/D  positions  of  2,  4,  6,  8,  10,  12,  16, 
and  20.  These  ratios  of  velocity  in  the  wake  to  velocity  in 
free  stream  were  extracted  from  pressure  coefficient  measure¬ 
ments  which  were  presented  in  Progress  Report  No  17  and  in  Ref  1 
These  same  velocity  distributions  are  presented  on  geometrically 
accurate  graphs  in  Pigs  1-2  through  1-13. 

For  purposes  of  this  analysis,  the  wake  boundary  is 
defined  by  the  points  where  the  velocity  defect  is  ten  per  cent 
of  the  velocity  defect  at  the  centerline,  i.e.,  (WVb>/(WvCL> 
=  0.10.  Using  this  definition,  experimental  wake  widths  were 
measured  for  each  X/D  position  for  the  twelve  bodies  above. 

These  experimental  wake  widths  are  indicated  in  Pigs  1-2  through 
1-13,  and  are  tabulated  in  Table  1-13. 

If  the  relationship  of  Eqn  (l.l)  is  to  be  an  accurate 
representation  of  the  wake  width,  it  is  necessary  that  the 
experimental  wake  widths  plot  as  straight  lines  on  log-log 
coordinate  paper.  Such  plots  are  presented  in  Pigs  1-14  through 
1-25  for  the  twelve  bodies  under  consideration.  It  is  seen  that 
for  the  more  or  less  streamlined  bodies,  in  particular.  Models 
P  through  L,  the  experimental  wake  widths  plot  as  a  straight 
line  within  the  limits  of  experimental  error.  This  is  also 
true  for  Model  A,  the  hollDw  hemisphere  open  to  the  flow.  For 
the  bluff  bodies.  Models  B  through  E,  the  experimental  wake 
widths  are  found  to  lie  on  one  straight  line  in  the  region 
2  ^  X/D  ^  8,  and  on  a  second  straight  line  in  the  region 
8  ^  X/D  N<  20. 

For  the  twelve  bodies  above,  an  approximate  represents 
tlon  of  wake  width  in  the  region  2  ^  X/D  ^  20  can  be  obtained 
by  assuming  a  straight  line  through  the  experimental  wake  widths 


3 


as  shown  in  Figs  1-14  through  l-25>  and  determining  the  empirical 
constants  of  Eqn  (l.l)  from  this  straight  line*  As  stated  above , 
the  accuracy  of  this  assumption  will  be  best  for  Model  A  and  the 
streamlined  bodies ,  Models  F  through  L,  although  Models  I  through 
L  show  some  large  deviations,  particularly  at  X/D  *  4.  It  is 
felt  that  this  deviation  is  caused  by  the  movement  of  the 
turbulent  boundary  layer  separation  point  downstream  to  some 
position  on  the  afterbody  of  these  models,  giving  an  initially 
narrower  wake.  Using  this  orte-line  approximation  of  the 
experimental  wake  widths,  the  empirical  constants  were  deter¬ 
mined  for  each  body,  and  are  presented  in  Table  1-14.  Using 
these  constants,  the  empirical  wake  widths  were  calculated  and 
are  shown  superimposed  on  the  velocity  distributions  in  Figs  1-2 
through  1-13.  It  is  seen  that  this  one  line  approximation  is 
indeed  very  accurate  for  Model  A  and  Models  F  through  L. 

Since  the  region  2  4  X/D  4  8  is  of  particular  interest 
in  problems  of  aerodynamic  retardation,  a  better  representation 
of  Wake  width  behind  each  of  the  bluff  bodies  (Models  B  through 
E)  can  be  obtained  by  assuming  one  straight  line  approximation 
in  the  region  2  4  X/D  4  8  and  a  second  in  the  region  8  ^  X/D  4  20. 
These  two-line  approximations  are  Bhown  in  Figs  1-15  through 
1-18,  and  the  constants  thus  determined  are  given  in  Table  1-14. 
The  empirical  wake  widths  then  calculated  are  superimposed  on 
the  lower  halves  of  Figs  1-3  through  1-6;  it  is  seen  that  such 
a  two-line  approximation  is  very  accurate  for  these  bodies. 

When  the  one -line  approximations  for  all  bodies  are 
plotted  together  as  in  Fig  1-26,  it  is  seen  that  three  general 
groupings  of  wake  width  according  to  the  body  shape  (or  Cp) 
are  Indicated.  As  shown  in  Fig  1-28  A,  wakes  close  behind 
bluff  bodies  similar  to  Models  A  through  E  can  be  represented 
in  a  first  approximation  by  a  single  straight  line  giving  the 
empirical  constants 


Figure  1-28  B  shows  the  single  approximating  line  for  inter¬ 
mediate  body  shapes  similar  to  Models  F  through  H,  with  the 
corresponding  wake  width  given  by 


b 

FT 


(l.lb) 


Similarly,  Fig  1-28  C  shows  the  single  approximating  line  for 
streamlined  bodies  similar  to  Models  I  through  L,  with  the 
corresponding  wake  width  given  by 


b 

IT 


(X.lc) 


It  is  seen  that  in  these  general  approximations  the  exponent  n 
is  nearly  the  same  for  the  three  body  groups.  The  effect  of 
increased  streamlining  of  the  bodies  or  of  a  decreasing  drag 
coefficient  is  a  decreasing  multiplier  k. 

When  the  two-line  approximations  shown  in  Figs  1-15 
through  1-18  are  plotted  together  as  in  Fig  1-27,  it  is  seen 
that  the  single  approximating  curve  shown  will  quite  accurately 
represent  the  wake  widths  of  Models  B  through  E.  Determining  the 
constants  for  this  general  approximation  yields  for  the  region 
2  4  X/D  ^  8  , 


b 


(l.ld) 


and  for  the  region  8  $  X/D  {  20  , 


b 


0.78 


X 

F 


0.63 


(l.le) 


2.1.3  Summary 


The  wake  widths  behind  twelve  bodies  of  revolution 
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In  subsonic  flow  have  been  measured  experimentally.  It  has  been 
found  possible  to  represent  these  wake  widths  for  all  bodies 
over  the  range  2  4  X/D  4  20  by  Eqn  (l.l),  with  the  appropriate 
constants  for  each  body  presented  in  Table  1-14.  For  the  bluff 
bodies  (Models  B  through  E)  a  better  representation  of  wake 
width  is  obtained  by  considering  separately  the  regions  2  4  X/D 4  8 
and  8  4  X/D  4  20;  the  appropriate  constants  are  also  given  in 
Table  1-14. 

A  more  general  representation  of  wake  width  is  obtained 
by  dividing  the  above  bodies  into  groups  of  bluff,  intermediate, 
and  streamlined  bodies,  with  Eqns  (l.la),  (l.lb),  and  (l.lc) 
giving  the  appropriate  empirical  constants.  These  general 
empirical  formulas  should  be  sufficiently  accurate  for  most 
wake  problems. 

2.1.4  Proposed  Work 

As  stated  In  Progress  Report  No  17,  Sec  2.1.5,  it  has 
not  been  possible  to  obtain  good  agreement  between  the  theoreti¬ 
cal  pressure  distributions  in  the  wake  as  predicted  by  Ref  2 
and  experimental  data.  In  the  region  close  behind  a  body.  The 
theory  uses  the  classical  assumption  that  the  wake  width,  b. 
Increases  with  the  1/3  power  of  X/D.  From  the  above  analysis 
it  is  seen  that  this  assumption  is  not  correct  In  the  region 
2  4  X/D  4  20  behind  bodies  of  revolution.  Therefore,  du^*ing 
the  next  reporting  period  the  theory  of  Ref  2  will  be  modified 
using  the  above  results  in  an  attempt  to  obtain  closer  agreement 
between  theoretical  and  experimental  pressure  distributions  In 
the  wake  close  behind  bodies  of  revolution. 

2.2  Theoretical  Analysis  of  Turbulent  Wake  In  Supersonic 

Flow 

2.2.1  Introduction 
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Progress  Report  No  19  presented  a  set  of  linearized 
partial  differential  equations  describing  the  velocity, 
pressure,  and  the  temperature  in  the  wake  of  an  axially  sym¬ 
metric  body.  During  the  last  reporting  period,  solutions  for 
the  set  of  equations  were  obtained  and  presented.  The 
solutions  are: 

u*  =  K1h(x*)e~h(x#)r#2, 

v*  =  2BK12f(x*)h3(x*)r#e'h^x*'r*2, 

p*  «  4B2K13fM<,2h4  (x*)f 2  (x*)  £  h(x*)r*2  + 

+  Ujxll  _  21  e-h(x*)r*2 

4BK1h2(x*)f2(x»)  J 


Where 


T«  =  Kt,h1(x*)e'hl^x#^r*2, 

r  f  1-1/2 

h(x*)  -  8BKi  I  f(x*)dx*  +  Kg  ' 
hl(x#)  “  [  +  K3  ]  ' 


with  terms  defined  as  in  Section  2.4,  Progress  Report  No  19. 

During  this  reporting  period  an  attempt  was  made  to 
determine  the  integration  constants,  Kj,  Kg,  K^,  and  K4.  The 
results  of  this  attempt  are  to  date  inconclusive. 

Since  the  solutions  of  the  equations  are  valid  for 
relatively  large  distances  downstream  of  the  body,  expanded 
experimental  data  for  the  larger  downstream  distances  is 
necessary  for  a  comparison  between  theory  and  experiment. 
Without  this  wider  range  of  data,  the  correlation  of  theory 
with  experiment  is  difficult  to  obtain. 
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Proposed  Work 


2.2.2 


Experimental  data  presently  available  for  shorter 
distances  downstream  of  various  bodies  will  be  used  In  further 
attempts  to  compare  theory  with  experimental  results.  In 
addition,  consideration  will  be  given  to  conducting  a  limited 
number  of  wake  study  experiments  at  larger  distances  downstream 
of  the  body  than  are  provided  for  by  present  data. 
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TABLE  VI.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  A 
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TABLE  1-2.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  B 
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1.00 

1.00 

0.99 

0.98 

1.00 

1.00 

1.00 

1.00 

0.99 

0.98 

0.95 

0.97 

1.00 

0.98 

0.99 

0.98 

0.97 

0.97 

0.96 

0.96 

2.00 

1.00 

0.95 

0.95 

0.95 

0.95 

0.95 

0.95 

0.96 

0.41 

0.82 

0.88 

0.91 

0.92 

0.94 

0.94 

0.95 

H  1 

0.14 

RH 

IamJ 

0.90 

0.92 

mm 

m 

0.95 

0.80 

0.72 

0.85 

0.90 

0.93 

0.94 

0.95 

m 

0.69 

0.85 

0.90 

0.90 

0.93 

0.94 

0.95 

mi 

0.91 

wm 

0.94 

0.95 

0 

0.66 

0.85 

0.89 

0.91 

0.92 

0.93 

0.94 

0.27 

0.67 

0.85 

0.69 

0.91 

0.93 

0.94 

0.95 

0.53 

0.69 

0.85 

0.90 

0.90 

0.93 

0.94 

0.95 

0.80 

0.72 

0.85 

0.90 

0.90 

0.93 

0.94 

0.95 

m 

0.14 

0.77 

0.87 

0.90  i 

0.92 

0.93 

0.94 

0.95 

1.33 

0.41 

0.82 

0.88 

0.91 

0.92 

0.94 

0.94 

0.95 

2.00 

1.00 

1  M 1 

0.95 

0.95 

0.95 

0.95 

0.95 

0.96 

2.67 

1.00 

0.98 

0.99 

0.98 

0.97 

0.97 

0.96 

0.96 

H  E  \  1 

1.00 

1.00 

1.00 

1.00 

0.99 

0.98 

0.95 

0.97 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.98 

4.67 

1.00 

1.00 

1.00 

1 . 00 

1.00 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

EKBI 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  ^.EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  C 


x/d 

2 

4 

6 

8 

10 

12 

16 

20 

fBl 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

3.33 

1.00 

0.99 

0.99 

1.00 

0.99 

0.99 

0.98 

0.98 

1.00 

0.99 

0.99 

0.99 

0.98 

0.98 

0.98 

0.98 

2.00 

1.00 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

1.33 

0.66 

0.85 

0.90 

0.92 

0.93 

0.94 

0.95 

0.96 

1.07 

0.41 

0.82 

0.89 

0.92 

0.92 

0.93 

0.94 

0.96 

v 

0.80 

0.17 

0.79 

0.88 

0.90 

0.92 

0.93 

0.94 

0.96 

y 

\4 

1RI 

0.17 

0.77 

0.87 

0.90 

0.92 

0.93 

0.94 

0.96 

0.27 

0.14 

0.75 

0.86 

0.90 

0.92 

0.93 

0.93 

0.96 

0 

0.14 

0.75 

0.86 

0.90 

0.92 

0.93 

0.93 

0.96 

0.27 

0.14 

0.75 

0.86 

0.90 

0.92 

0.93 

0.93 

0.96 

0.53 

0.17 

0.77 

0.87 

0.90 

0.92 

0.93 

0.94 

0.96 

0.80 

0.17 

0.79 

0.88 

0.90 

0.92 

0.93 

0.9^ 

0.96 

0.41 

0.82 

0.89 

0.92 

0.92 

0.93 

0.94 

0.96 

0.66 

0.85 

0.90 

0.92 

0.93 

0.94 

0.95 

0.96 

2.00 

1.00 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

0.96 

1.00 

0.99 

0.99 

0.98 

0.98 

0.98 

0.98 

0.98 

3.33 

1.00 

0.99 

0.99 

1.00 

0.99 

0.99 

0.98 

0.98 

4.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

am 

1.00 

!  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

mu 

BiH 

IBI 1 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-4.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  D 


>l>' 


x/d 

2 

4 

6 

8 

10 

12 

16 

20 

v 

7.67 

1.00 

1.00 

1.00 

1  1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1  1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1.00 

;  1.00 

1.00 

1.00 

1.00 

IHH 

4.00 

1  1.00 

1.00 

1.00 

0.99 

0.99 

rai 

mm 

3.33 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.99 

0.98 

2.67 

0.94 

0.9S 

0.99 

0.98 

0.98 

0.98 

0.97 

0.97 

2.00 

0.94 

0.97 

0.97 

0.96 

0.96 

0.96 

0.95 

0.96 

1.33 

0.68 

0.&7 

0.90 

0.92 

0.93 

:  0.94 

0.95 

0.95 

1.07 

0.40 

0.82 

0.88 

0.90 

0.92 

ran 

rai 

0.80 

0.17 

0.77 

0.87 

!  0.90 

0.92 

raa 

raa 

V 

V. 

0.53 

0.00 

0.73 

0.87 

0.90 

0.92 

0.94 

0.95 

0.95 

0.27 

0.00 

0.71 

0.86 

0.90 

0.92 

0.94 

0.95 

0.95 

0 

0.00 

0.71 

0.86 

0.90 

0.91 

0.93 

0.94 

0.95 

0.27 

j  0.00 

0.71 

0.86 

0.90 

0.92 

0.94 

0.95 

0.95 

0.53 

;  0.00 

0.73 

0.87 

0.90  1 

0.92 

0.94 

0.95 

0.95 

0.80 

0.17 

0.77 

0.87 

0.90 

0.92 

0.94 

0.95 

0.95 

1.07 

0.40 

0.82 

0.88 

0.90  ; 

0.92 

0.94 

0.94 

0.95 

1.33 

0.68 

0.87 

0.90 

0.92 

0.93 

0.94 

0.95 

0.95 

2.00 

0.94 

0.97 

0.97 

0.96 

0.96 

0.96  i 

0.96 

0.96 

2.67 

0.94 

0.99 

0.99 

0.98 

0.98 

0.98 

0.97 

0.97 

3.33 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.98 

O.98 

4.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.99 

O.98 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  . 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-5. EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  E 


5* 

vb 

2 

4 

6 

8 

10 

12 

16 

20 

7.67 

1.00 

1.00 

H* 

• 

0 

0 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

wm 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

M 

• 

O 

O 

I  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

IRS 

1.00 

n 

1.00 

1.00 

1.00 

1 1 

2.67 

• 

0 

0 

1.00 

1.00 

1.00 

0.99 

0.99 

0.98 

0.98 

2.00 

0 

0 

• 

rH 

0.99 

0.98 

0.97 

0.97 

0.97 

0.97 

0.97 

1.33 

0.92 

0.93 

0.94 

0.94 

0.94 

1  9 :  9B 

0.96 

0.96 

1.07 

0.73 

0.87 

0.91 

0.93 

0.94 

0.94 

0.96 

0.96 

.<l< 

0.80 

0.47 

0.84 

0.90 

0.92 

0.93 

0.94 

0.96 

0.96 

0.53 

0.20 

0.81 

0.90 

0.92 

0.93 

0.94 

0.96 

0.96 

0.27 

0.00 

0.80 

0.88 

0.92 

0.93 

0.94 

0.95 

0.96 

0 

0.00 

0.80 

0.88 

0.91 

IH9I 

0.94 

0.95  1 

0.96 

0.00 

0.88 

0.92 

m 

0.94 

wm 

0.96 

0.53 

0.20 

0.81 

0.90 

0.92 

0.93 

0.94 

0.96 

0.96 

0.80 

0.47 

0.84 

0.90 

0.92 

0.93 

0.94 

0.96  i 

0.96 

1.07 

0.73 

0.67 

0.91 

0.93 

0.94 

0.94 

0.96 

0.96 

1.33 

0.92 

0.91 

0.94 

0.94 

0.94 

0.94 

0.96 

0.96 

2.00 

1.00 

0.99 

0.98 

0.97 

0.97 

0.97 

0.97 

0.97 

2.67 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.98 

0.98 

3.33 

1.00 

1.00 

1.00 

1.00 

wm 

1.00 

wm 

■99 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

'1.00 

1.0c 

7.67 1 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-6.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  F 


x/fe 


7.67 


5.67 


4.67 


4.00 


3.33 


2.67 


2.00 


1.33 


1.07 


0.60 


0.53 


0.27 


0 


0.27 


0.53 


0.60 


1.07 


1.33 


2.00 


2.67 


.33 


4.00 


4.67 


5.67 


7.67 


2 

4 

6 

8 

10 

12 

16 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

;  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.96 

1.00 

0.99 

0.99 

1.00 

1.00 

0.96 

0.95 

0.95 

0.96 

0.98 

0.58 

0.93 

0.95 

0.95 

0.95 

0.94 

0.95 

0.96 

0.97 

0.82 

0.90 

0.93 

0.94 

0.94 

0.95 

0.96 

0.96 

0.68 

0.88 

0.50 

0.93 

0.93 

0.95 

0.95 

0.96 

0.50 

0.82 

0.88 

0.92 

0.93 

0.95 

0.95 

0.96 

0.33 

0.81 

0.88 

0.90 

0.92 

0.94 

0.95 

0.95 

0.29 

0.80 

0.87 

0.90 

0.93 

0.93 

0.94 

0.95 

0.33 

0.81 

0.88 

0.90 

0.93 

0.94 

0.95 

0.95 

0.50 

0.82 

0.88 

0.92 

0.94 

0.95 

0.95 

0.96 

0.68 

0.88 

0.90 

0.93 

0.94 

0.95 

0.95 

0.96 

0.82 

0.90 

0.93 

0.94 

0.95 

0.95 

O.96 

0.96 

0.93 

0.95 

0.95 

0.95 

0.9£ 

0.95 

0.56 

0.97 

1.00 

1.00 

0.98 

0.99 

0.99 

0.98 

0.96 

0.98 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1 . 00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.0c 

1.00 

1.00 

1.0 

1.0c 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-7  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  G 


^R 

X/D 

2 

4 

6 

8 

10 

12 

16 

20 

7.67 

1.00 

1.00 

s  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1  1.00 

!  1.00 

1.00 

1.00 

0.99 

0.99 

4.00 

0.99 

1.00 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

3.33 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0 

VO 

0 0 

2.67 

0.99 

0.99 

0.99 

0.99 

0.99 

O.98 

0.98 

0.96 

2.00 

1HHI 

0.99 

0.99 

0.99 

0.98 

O.96 

0.98 

0.97 

m 

0.96 

0.98 

0.98 

0.97 

0.97 

0.97 

0.98 

0.97 

0.97 

0.97 

0.97 

0.96 

O.96 

0.96 

0.97 

0.96 

0.96 

HH1 

0.95 

0.95 

0.95 

O.96 

0.96 

0.96 

v« 

0.53 

0.89 

0.92 

0.93 

0.95 

0.94 

0.94 

0.95 

0.95 

0.27 

0.79 

0.90 

0.91 

0.92 

0.93 

0.94 

0.95 

0.95 

0 

0.74 

0.88 

0.90 

0.93 

0.93 

0.94 

0.95 

0.95 

0.90 

0.91 

0.92 

0.93 

0.94 

0.95 

0.95 

m 

0.^9 

0.92 

0.93 

0.95 

0.94 

0.94 

0.95 

0.95 

0.60 

0.96 

0.95 

0.95 

0.95 

0.95 

O.96 

0.96 

0.96 

HR 

EBB 

0.97 

0.97 

0.96 

0.96 

O.96 

0.97 

0.96 

H9 

0.96 

0.98 

0.97 

0.97 

0.97 

0.9& 

0.97 

2.00 

0.99 

0.99 

0.99 

0.99 

0.98 

0.98 

0.96 

0.97 

2.67 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

0.98 

0.98 

3.33 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

4.00 

0.99 

1.00 

0.99 

0.99 

0.99 

0.99 

0.99 

0.98 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

HS 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-8.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  H 


r/R 

Vu 

2 

4 

6 

8 

10 

12 

16 

20 

<<l< 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  ’ 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.99 

0.99 

1.33 

0.99 

0.99 

0.97 

0.96 

0.97 

0.96 

0.96 

0.96 

1.07 

0.98 

0.97 

0.95 

0.95 

0.95 

0.96 

0.95 

0.96 

0.80 

0.93 

0.93 

0.93 

0.94 

0.94 

0.94 

0.95 

0.96 

0.53 

0.84 

0*88 

0.90 

0.92 

0.93 

0.93 

0.94 

0.96 

0.27 

0.75 

0.84 

0.88 

0.92 

0.92 

0.95 

0 

0.74 

0.81 

0.86 

0.90 

0.90 

0.92 

0.93 

0.95 

0.27 

0.75 

0.84 

0.88 

0.92 

0.92 

0.92 

0.93 

0.95 

0.53 

0.84 

0.88 

0.90 

0.92 

0.93 

0.93 

0.94 

0.96 

0.80 

0.93 

0.93 

0.93 

0.94 

0.94 

0.94 

0.95 

0.96 

1.07 

0.98 

0.97 

0.95 

0.95 

0.95 

0.96 

0.95 

0.96 

1.33 

0.99 

0.99 

0.97 

0.96 

0.97 

0.96 

0.96 

0.96 

2.00 

1.00 

1.00 

1.00 

1.00 

Ui'v3:] 

•  '  1 

■M/LJ 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

;  1.00 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-9.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  I 


2 

4 

6 

8 

10 

12 

16 

20 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  ! 

1.00 

4.67 

:.oo 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

2.00 

1.00 

1.00 

o 

o 

• 

1.00 

1.00 

1.00 

0.99 

0.98 

1.33 

o.99 

1.00 

0.97 

0.97 

0.97 

0.97 

0.96 

0.96 

1.07 

0.99 

0.99 

0.94 

0.95 

0.95 

0.95 

0.95 

0.95 

V 

0.80 

0.98 

0.95 

0.90 

0.93 

0.93 

0.93 

0.95 

0.96 

V. 

0.53 

0.84 

i  0.86 

0.90 

0.90 

0.91 

0.91 

0.94 

0.95 

liH 

0.71 

0.79 

0.89 

0.86 

0.90 

0.90 

0.93 

0.95 

0 

0.66 

0.75 

0.88 

0.85 

0.89 

0.90 

0.92 

0.94 

0.27 

0.71 

0.79 

0.89 

0.86 

0.90 

0.90 

0.93 

0.95 

0.53 

0.84 

0.86 

0.90 

0.90 

0.91 

0.91 

0.94 

0.95 

0.80 

0.98 

0.95 

0.90 

0.93 

0.93 

0.93 

0.95 

0.96 

1.07 

0.99 

0.99 

0.94 

0.95 

0.95 

0.95 

0.95 

0.95 

1.33 

0.99 

1.00 

0.97 

0.97 

0.97 

0.97 

0.96 

O.96 

2.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

O.98 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1^1 

0.99 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.0C 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

e m 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  1-10.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  J 


>l>* 


%  I _ x/p 


2 

4 

6 

8 

10 

12 

16 

20 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

S.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

;  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

l.C 

4.00 

1.00 

1.00 

1  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

:  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2.67 

1.00 

!  1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

2.00 

0.99 

:  1.00 

1  1.00 

1.00 

0.99 

1.00 

0.99 

0.99 

1.33 

0.99 

0.99 

0.99 

0.98 

0.98 

0.97 

0.96 

0.97 

1.07 

0.99 

0.99 

0.97 

0.95 

mm 

0.80 

0.98 

0.96 

0.93 

1  0.92 

H 

wm 

o.9» 

0.94 

0.53 

0.87 

0.85 

0.68 

0.91 

0.90 

0.92 

0.94 

0.27 

0.72 

0.76 

0.83 

0.86 

0.86 

0.90 

0.93 

0.93 

0 

0.64 

0.71 

0.83 

0.86 

0.85 

0.68 

0.92 

0.92 

0.27 

0.72 

0.76 

O.83 

0.86 

0.86 

0.90 

0.93 

0.93 

0.S3 

0.87 

0.85 

0.88 

0.91 

0.90 

0.92 

0.93 

0.94 

0.80 

O.98 

0.96 

0.93 

0.92 

0.93 

0.93 

0.94 

0.94 

1.07  | 

0.99 

0.99 

0.97 

0.95 

0.96 

0.95 

0.95 

0.95 

1.33 

0.99 

0.99 

0.99 

0.98 

0.98 

0.97 

0.96 

0.97 

2.00 

bbihbi  1 

1.00 

1.00 

1.00 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00  , 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  : 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  ! 

1.00 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

TABLE  Ml.  EXPERIMENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  K 
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2 

4 

6 

8 

10 

7.67 

1.00 

1.00 

1.00 

1.00 

1.00 

5.67 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

;  1.00 

1.00 

1.00 

1.00 

1.00 

4.00  S 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

!  1.00 

1.00 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

2.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.33 

0.95 

0.99 

0.98 

0.97 

0.98 

1.07 

0.99 

0.99 

0.95 

0.94 

0.96 

0.80 

0.97 

0.96 

0.90 

0.91 

0.93 

0.53 

0.68 

0.90 

0.90 

0.89 

0.90 

0.27 

0.72 

0.77 

0.87 

0.88 

0.86 

'  0 

0.66 

0.71 

0.88 

0.89 

0.83 

0.27 

0.72 

0.77 

0.87 

0.88 

0.86 

0.53 

0.88 

0.90 

0.90 

O.89 

0.90 

0.80 

0.97 

O.96 

0.90 

0.91 

0.93 

0.99 

0.99 

0.95 

0.94 

0.96 

1^9 

mm 

mm 

R 1  - 

.  |  "'‘"I 

1.00 

1.00 

1.00 

1.00 

2.67 

1.00 

1.00 

1.00 

1.00 

1.00 

3.33 

1.00 

1.00 

1.00 

1.00 

1.00 

4.00 

1.00 

1.00 

1.00 

1.00 

1.00 

4.67 

1.00 

1.00 

1.00 

1.00 

1.00 

Rfl 

1.00 

1.00 

1.00 

1.00 

1.00 

7.67 1 

1.00 

1.00 

1.00 

1.00 

1.00 

12 

16 

20 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

0.99 

0.99 

0.99 

0.98 

0.95 

0.96 

0.95 

0.94 

0.95 

0.93 

0.94 

0.95 

0.91 

0.94 

0.94 

0.90 

0.93 

0.93 

0.90 

0.93 

0.93 

0.90 

0.93 

0.93 

0.91 

0.94 

0.94 

0.93 

0.94 

0.95 

0.95 

0.94 

0.95 

0.98 

0.96 

0.99 

0.99 

0.99 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00  i 

1.00  j 

1.00 

1.00 

1.00 

1.00 

TABLE  T12.EX  PERI  MENTAL  VELOCITY  DISTRIBUTION 

IN  WAKE  OF  MODEL  L 


MODEL 


*/d=2 


12  16 


b/R 

1.70 

2.55 

3.50 

3.70 

4.00 

4.50 

6.00 

6.30 

1.88 

2.20 

2.60 

3.05 

3.70 

4.10 

4.60 

5.40 

*>t» 

1.75 

2.25 

2.55 

2.90 

3.40 

4.00 

4.50 

1.70 

2.30 

2.60 

2.80 

3.40 

0 

CM 

• 

4.80 

*Vr 

1.70 

2.15 

2.50 

2.80 

3.10 

3.80 

4.60 

5.00 

^■r 

1.30 

1.85 

2.25 

2.50 

2.80 

2.90 

3.50 

0 

CM 

• 

^■r 

1.35 

1.75 

2.25 

2.25 

2.70 

2.90 

2.90 

3.40 

t?/R 

1.10 

1.70 

2.20 

3.00 

3.30 

3.60 

4.30 

5.00 

Vr 

1.03 

1.20 

B5 1 

1.95 

2.00 

2.00 

2.30 

2.50 

b/R 

0.80 

0.98 

1.65 

1.55 

1.75 

lali 

2.40 

2.80 

b/R 

0.72 

0.85 

1.20 

1.40 

1.60 

1.90 

2.20 

2.30 

b/fc 

0.82 

0.90 

1.50 

1.75 

1.70 

1.70 

2.20 

2.50 

TABLE  1-13.  EXPERIMENTAL  WAKE  WIDTHS 


MODEL 


c.. 


1.40 


2<X/D<20|  2<X/D<8  |  8<  X/D<20 


1.14  0.582 


!  1.11 

1.26 

0.454 

1.39 

0.368 

0.788 

0.642 

1.03 

1.33 

0.394 

1.36 

0.365 

0.734 

0.660 

O.836 

1.34 

0.384 

tm 

0.879 

0.574 

0.743 

1.32 

0.393 

1.33 

0.358 

0.684 

0.678 

0.577 


0.425 


0.350 


0.230 


0.206 


0.197 


0.194 


0.488 


O.665 


o.8iq  0.376 


0.522 


0.531 


H 


0.565  0.491 


TABLE  1-14.  EMPIRICAL  WAKE  WIDTH  CONSTANTS 
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FIG  1-10. EXPERIMENTAL  VELOCITY  DISTRIBUTION  IN  WVKE  OF  OGIVE  CYLINDER  WITH  AFTERBODY  NO.  1 - MODEL 
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FIG  FI 2  EXPERIMENTAL  VELOCITY  DISTRIBUTION  IN  WAKE  OF  OGIVE  CYLINDER  WITH  AFTERBODY  NQ  3 — MODEL 
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3.0  Investigation  of  Basic  Stability  Parameters  of 

Conventional  Parachutes 

3.1  Experimental  Determination  of  the  Apparent  Moment 
of  Inertia  of  Parachutes 

3.1.1  Introduction 

In  Progress  Report  No  20,  an  experimental  method  to 
determine  the  apparent  moment  of  inertia  was  presented.  The 
theory  and  supporting  references  were  also  presented  in  the 
same  report. 

During  the  present  reporting  period,  tests  were  con¬ 
ducted  with  simple  geometric  bodies.  The  results  obtained  were 
compared  with  those  of  the  references  to  check  the  validity  of 
the  experimental  procedure  used  and  the  accuracy  of  the  test 
results . 

3.1.2  Test  Apparatus  Modifications 

It  was  found  during  the  initial  tests  that  the 
mounting  of  the  "Shaevitz"  transducer  was  not  rigid  enough. 
Since  the  allowable  tolerance  in  the  setting  of  the  core  and 
the  coil  of  the  transducer  is  very  small,  the  core  was  attached 
to  a  more  rigid  mounting.  With  this  modification,  the  system 
performs  satisfactorily. 

By  attaching  a  pulley  to  the  test  frame,  it  was 
possible  to  use  static  loading  (Pig  4-l)  as  a  means  of  deter¬ 
mining  the  spring  constant. 
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3-1-3 


Models 


In  the  preliminary  experiments  conducted  during  this 
reporting  period,  the  models  used  were  similar  to  those  used 
In  Refs  1  and  2.  These  models  were  used  to  enable  comparison 
of  recent  experimental  results  obtained  with  the  modified 
apparatus  with  the  values  obtained  by  previous  experimenters 
and  by  theoretical  methods.  The  three  bodies  used  in  the  experi¬ 
ments  consisted  of  a  sphere,  a  cube,  and  a  circular  disk. 
Reference  2  provides  experimental  data  for  cubes  and  spheres 
and  cites  further  references  for  theoretical  calculations  deal¬ 
ing  with  spheres .  Reference  1  contains  experimental  data  for 
several  disks  and  cites  theoretical  references  for  the  purpose 
of  comparing  theoretical  and  experimental  results. 

It  was  discovered  in  earlier  tests  that  the  size  of 
the  models  had  to  be  relatively  large  so  as  to  reduce  the 
relative  effects  of  the  system  friction.  The  optimum  sized 
model  for  the  system  used  has  a  mass  comparable  to  or  greater 
than  a  one  inch  diameter  brass  sphere. 

3.1.4  Experimental  Procedure 

The  test  models  are  mounted  on  the  apparatus  as  shown 
in  Fig  4-2  to  insure  symmetry  and  to  optimize  the  amplitude  of 
oscillation. 

Progress  Report  No  20,  Sec  3.1. 3*  presented  a  detailed 
description  of  the  basic  apparatus  used.  During  testing  of  the 
models,  the  oscillator  frequency  is  constantly  checked  for 
accuracy  by  feeding  the  signals  into  a  y-axJs  amplifier  of  an 
oscilloscope  and  feeding  60  cycle  line  current  into  the 
x-ampllfier.  By  observing  the  Lissajous  figure  appearing  on 
the  oscilloscope  screen,  any  whole  number  or  simple  fraction 
multiple  of  60  cycles  per  second  can  be  set  on  the  oscilloscope. 


In  these  tests,  a  frequency  of  15  cycles  per  second  was  used  as 
the  time  base,  corresponding  to  a  4  to  1  Lissajous  figure  on 
the  screen.  This  time  base  frequency  is  fed  into  one  channel 
of  a  Brush  Recorder.  Into  the  other  channel  is  fed  the  output 
of  the  strain  analyzer  which  in  turn  is  excited  by  the  trans¬ 
ducer  mounted  on  the  torsion  rod  system. 

The  determination  of  the  spring  constant  can  be  accom¬ 
plished  by  means  of  a  static  or  a  dynamic  test.  In  the  static 
test,  the  system  was  set  up  as  shown  in  Pig  4-1,  with  the  angle 
of  deflection  determined  by  the  small  angle  formula 

e-AY/<  ,  (n.i) 

Since  the  moment  arm  used  in  this  experiment  seemed 
to  bend  appreciably,  thereby  adding  its  own  deflection,  another 
test  was  conducted  to  correct  for  the  additional  deflection  of 
the  moment  arm.  The  moment  arm  was  secured  in  a  rigid  vice  and 
reloaded  in  the  same  manner  as  before.  The  deflection  for  any 
given  load  measured  in  this  experiment  was  subtracted  from  the 
deflection  for  the  corresponding  load  on  the  torsion  system. 

The  angles  were  plotted  as  a  function  of  the  applied  torque, 
with  the  slope  of  the  line  determining  the  spring  constant. 

In  the  dynamic  test,  the  moment  of  inertia  of  the 
system  was  calculated  and  the  period  in  air  measured.  Prom 
these  values  the  spring  constant  was  calculated  using  Eqn  4.2 
of  Progress  Report  No  20} 

K  sr  4ft2  ^  (4.2) 

?d 

where  K  is  the  spring  constant. 

The  calculated  moment  of  inertia  included  Systems 
1  and  2  shown  in  Pig  4-2.  System  3  gave  negligible  values  as 
compared  to  Systems  1  and  2.  System  1,  which  included  the 
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models  and  the  members  used  to  attach  the  models  to  the  torsion 

rod,  had  a  calculated  moment  of  inertia  of  60.8  x  10-5  gm-cm  . 

System  2,  which  included  the  fastener  and  the  rods  supporting 

the  detection  core,  had  a  calculated  moment  of  inertia  of 
q  2 

0.331  x  10  gm-cm  .  The  torsion  rod  itself,  which  had  a  moment 

o  p 

of  inertia  of  0.00049  x  10-3  gm-cm  ,  constituted  System  3. 

The  spring  constant  determined  in  this  manner  compared 
very  well  with  the  constant  determined  by  the  static  test.  The 
average  value  of  the  spring  constant  obtained  by  the  static  test 
is  45.8  x  10^  dyne-cm.  The  spring  constant  determined  by 
dynamic  tests  is  being  used  in  this  report  because  the  testing 
procedure  used  for  determining  the  added  mass  is  similar  to  the 
dynamic  determination  of  the  spring  constant. 

A  sample  of  recorder  paper  from  tests  on  the  1.5  inch 
diameter  sphere  is  shown  in  Fig  4-3.  The  figure  illustrates 
the  procedure  for  calculating  the  period  of  oscillation  of  the 
system. 

3.1.5  Test  Results 

The  data  from  tests  on  the  disk  is  presented  in  terms 
of  the  additional  mass,  M1,  in  order  that  results  may  be  com¬ 
pared  with  those  of  Ref  1.  The  working  equation  used  for  cal¬ 
culating  M'  was  arrived  at  by  modification  of  Eqn  (4.6)  of 
Progress  Report  No  20.  The  modified  equation  has  the  form 

I1  =  (Tw£  -  TA2)  -  IB.  (4.3) 

where  the  subscripts  w  and  A  refer  to  water  and  air  respectively. 

Under  present  testing  conditions,  the  added  moment  of 
inertia  of  the  system  without  models,  Is,  was  found  to  be  negli¬ 
gible.  Also,  for  the  whole  system t 

I«  «=  M«|^  .  (4.4) 
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Then  for  one  model , 


M« 


(4.5) 


To  facilitate  the  comparison  of  test  results  with  the 
results  of  Ref  2,  the  data  from  tests  on  the  cube  and  sphere 
models  Is  presented  as  the  dimensionless  ratio  where 

Is  the  mass  of  the  water  displaced  by  the  model.  These 
experimental  results  are  presented  In  Table  4-1,  where  column  1 
gives  the  characteristic  dimension  of  the  model  in  metric  units. 
Figure  4-4  is  a  graph  comparing  data  from  tests  on  disks  with 
experimental  and  theoretical  values  obtained  from  Ref  1.  This 
project's  test  results  lie  between  previously  obtained  experi¬ 
mental  values  and  the  values  calculated  from  available  theory. 
The  comparison  of  this  project's  experimental  results  and  Lamb’s 
theoretical  values  for  spheres  and  cubes  is  (presented  in 
Fig  4-5. 

Since  the  project's  test  results  compare  well  with 
previous  results  and  available  theory,  it  is  felt  that  the 
experimental  apparatus  used  and  the  test  procedure  developed 
are  suitable  for  the  experimental  determination  of  the  apparent 
mass  and  apparent  moment  of  inertia  of  rigid  canopy  models . 


3.1.6  Proposed  Work 


During  the  next  reporting  period,  models  of  circular 
flat,  ribless  guide  surface,  and  ribbon  parachute  canopies  will 
be  constructed  and  tested.  Further  refinements  in  the  experi¬ 
mental  apparatus  and  technique  will  be  made  where  necessary. 


3.2  Wind  Tunnel  Investigation  of  Parachute  Models 

A  second  draft  of  the  technical  report  entitled 
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"Stablillty  and  Drag  of  Parachutes  with  Varying  Effective 
Porosity,"  has  been  submitted  to  the  Procuring  Agency  for 
approval . 

3.3  Theoretical  Investigation  of  the  Dynamic  Stability 

of  Parachutes 

A  preliminary  technical  report  of  this  data  has  been 
written  and  is  currently  being  reviewed. 
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TABLE  4-1.  DATA  AND  CALCULATIONS  FOR  m'  AND  d  OF  VARIOUS  MODELS 
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4,0  Theoretical  Study  of  Supersonic  Parachute  Phenomena 

4.1  Introduction 

This  study  has  analyzed  the  behavior  of  conventional 
parachutes  in  supersonic  flow  in  the  wind  tunnel  and  by  means 
of  the  water  analogy.  Based  on  this  analysis,  a  supersonic 
decelerator  was  conceived  which  utilizes  a  pressure  reducing 
"spike"  in  connection  with  a  solid  cloth  canopy. 

As  originally  noted  in  Progress  Report  No  20,  the 
investigation  of  the  spike  concept  in  combination  with  a  ribbon 
parachute  canopy  will  be  presented  under  this  study.  It  was 
previously  presented  as  Project  No  10  of  Progress  Reports  Nos 
18  and  19. 

The  work  done  during  this  reporting  period  will  be 
reported  in  the  following  phases: 

a)  Supersonic  wind  tunnel  studies  of  flexible  spiked 
parachutes 

b)  Pressure  distribution  studies  on  spiked  parachutes 
in  supersonic  flow 

c)  Supersonic  wind  tunnel  studies  of  spiked  ribbon 
parachutes . 

4.2  Supersonic  Wind  Tunnel  Studies  of  Flexible  Spiked 
Parachutes 

4.2.1  Past  Work 

Rigid  models  of  the  spiked  solid  canopy  parachute, 
both  with  and  without  suspension  lines,  have  been  found 
stable  in  wind  tunnel  tests  at  Mach  numbers  of  1.14,  2.0,  and 
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3.0.  Test  results  for  configurations  with  textile  canopies 
were  presented  in  Progress  Reports  Nos  17 ,  18,  19,  and  20. 

Mechanical  damage  to  the  model  and  to  the  wind 
tunnel  suspension  system  occurring  in  the  early  tests  of  flex¬ 
ible  canopy  models  prompted  the  improvement  of  the  suspension 
line  connection,  design  of  a  new  deployment  device,  and  minor 
improvements  in  the  parachute  model.  The  new  deployment  device 
appeared  in  Pig  7-2  of  Progress  Report  No  20.  Damage  to  the 
model  was  reduced  through  the  use  of  heavier  MIL-C-8021A,  Type 
II  nylon  cloth  (300  lb/in),  and  by  elimination  of  the  wire  ring 
used  in  the  canopy  edges.  In  addition,  cones  are  now  molded  of 
low  density  "Arothane,"  a  product  of  the  Archer-Daniels -Midland 
Company  of  Minneapolis,  Minnesota.  These  lighter  weight  cones 
have  reduced  the  number  of  suspension  line  failures  while  at 
the  same  time  improving  the  dynamic  balance  of  the  system  as 
reported  in  Progress  Report  No  20. 

4.2.2  Continuation 

Tests  have  been  continued  during  the  present  reporting 
period  on  flexible  models  of  the  spiked  solid  canopy  parachute. 
The  aim  of  these  tests  is  to  further  improve  structural  reli¬ 
ability  and  strength,  and  to  establish  an  optimum  geometry  for 
the  configuration.  Improvements  made  in  the  suspension  lines, 
models,  and  deployment  device  continue  to  serve  very  satisfac¬ 
torily.  Figure  7-1  shows  a  detailed  canopy  edge  view  of  how 
suspension  lines  are  fastened  to  the  canopy  fabric. 

Continued  good  test  results  indicate  that  the  present 
spiked  solid  cloth  canopy  parachute  model  Is  very  near  to  an 
optimum  geometric  form.  During  the  present  reporting  period 
85#  of  the  wind  tunnel  tests  conducted  were  completely  success¬ 
ful.  These  results  and  the  valuable  Information  gained  from 


51 


the  limited  success  of  the  remaining  15$  of  the  tests  firmly 
indicate  the  good  repeatability  of  the  parachute's  performance. 

A  Schlleren  photograph  of  the  flexible  model  during  a  wind 
tunnel  test  is  shown  in  Fig  7-2. 

Test  measurements  to  determine  the  drag  of  the  flexible 
parachute  models  were  obtained  during  this  reporting  period. 

The  drag  coefficient,  based  on  total  cloth  area,  was  Cp  -  0.335. 
While  the  results  from  such  a  very  limited  investigation  are 
far  from  concilia ive,  the  drag  coefficient  obtained  is  in  the 
general  range  of  the  values  characteristic  of  supersonic  decel¬ 
eration  devices.  These  tests  incorporated  the  drag  balance 
arrangement  shown  in  Fig  7-3* 

In  the  latter  part  of  this  reporting  period,  two  wind 
tunnel  tests  were  conducted  on  a  flexible  parachute  model  in 
combination  with  an  ogive  forebody.  The  ogive  has  a  calibre 
of  2.5,  a  fineness  ratio  of  2.5*  and  a  base  diameter  of 
1.83  inches.  Size  and  distance  relationships  between  the  fore¬ 
body  and  the  parachute  model  are  Di/fl>2  •  0.458  and  «  9.75» 
where  D2  is  the  maximum  diameter  of  the  canopy,  the  diameter 
of  the  forebody,  and  L  is  the  distance  between  the  leading  edge 
of  the  canopy  and  the  base  of  the  forebody.  The  test  arrange¬ 
ment  used  in  these  wind  tunnel  tests  is  shown  in  Fig  7-4. 

Further  testing  will  be  necessary  before  any  definite  results 

concerning  the  forebody  experiments  can  be  presented. 

% 

4.2,3  Proposed  Work 

Supersonic  wind  tunnel  tests  will  be  continued  on  the 
spiked  sc.  id  cloth  canopy  models  to  insure  an  optimum  model 
geometry.  Drag  measurements  in  free  stream  will  be  continued 
in  order  to  obtain  more  conclusive  results.  In  addition,  test¬ 
ing  will  be  conducted  to  determine  the  effect  of  the  forebody 
on  the  drag  of  the  flexible  spiked  parachute. 
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4.3  Pressure  Distribution  Studies  on  Spiked  Parachutes 
in  Supersonic  Plow 

Progress  Report  No  20  presented  tabulated  results  from 
the  pressure  distribution  studies  of  rigid  models  of  the  spiked 
solid  canopy  parachute.  Table  7-1  and  Pig  7-4  of  that  report 
presented  the  data  from  Mach  3  tests  In  terms  of  pressure 
coefficients  as  It  Is  defined  In  Sec  4.3.3  of  Progress  Report 
No  19.  No  further  work  In  this  direction  has  been  undertaken 
during  the  present  reporting  period. 

It  is  proposed  that  during  the  next  reporting  period 
pressure  distribution  studies  on  rigid  models  of  the  spiked 
solid  cloth  canopy  parachute  in  the  wake  of  a  forebody  be  con¬ 
ducted.  This  testing,  in  conjunction  with  the  forebody  studies 
on  flexible  models  of  the  spiked  parachute  as  noted  previously, 
will  make  it  possible  to  obtain  a  parachute  design  that  performs 
satisfactorily  in  the  wake  of  a  forebody. 

In  conjunction  with  the  previously  described  test 
program,  a  series  of  tests  are  also  being  conducted  at  the  Naval 
Ordnance  Laboratory  at  Silver  Springs,  Maryland.  This  particular 
segment  of  the  test  program  is  under  the  direction  of  Mr.  George 
Sloan.  In  these  tests,  a  two  Inch  diameter  steel  model  of  the 
spiked  solid  canopy  parachute  is  towed  behind  a  heavy  forebody 
which  is  fired  from  a  cold  gas  gun.  A  description  of  the  test 
facilities  in  Aerodynamics  Range  No  1  used  for  these  tests  can 
be  found  In  the  Naval  Ordnance  laboratory  publication, 
AEROBALLISTIC  RESEARCH  FACILITIES,  NOLR  1233.  At  the  end  of 
this  reporting  period  only  preliminary  tests  had  been  aonducted. 

V 

Figure  7-5  shows  a  typical  model  and  forebody  with  connecting 
cable. 

4.4  Supersonic  Wind  Tunnel  Studies  of  Spiked  Ribbon 
Parachutes 


4*4*1 


Introduction 


The  "spike”  concept  used  successfully  In  the  previously 
mentioned  parachute  has  been  used  In  conjunction  with  ribbon 
parachute  canopies.  The  water  analogy  tests  conducted  on  this 
type  of  configuration,  as  presented  In  Progress  Report  No  18, 
Project  10,  lead  to  the  construction  of  rigid  models  o t  the 
parachute.  Results  from  tests  on  rigid  models  with  geonetric 
porosities  of  20,  25,  and  30  per  cent  distributed  evenly  over 
the  canopy  were  originally  presented  In  Table  7-2  of  Progress 
Report  No  20.  During  this  reporting  period,  tests  have  been 
conducted  on  a  rigid  model  with  40#  geometric  porosity  distri¬ 
buted  evenly  over  the  canopy  (Pig  7-6). 

4.4.2  Tests 

The  rigid  models  are  tested  In  a  supersonic  wind  tunnel 
at  Mach  3.0.  Figure  7-7  of  Progress  Report  No  20  shows  the  test 
arrangement  used.  During  the  test  runs,  visual  observation  and 
high  speed  Schlieren  movies  are  used  to  determine  the  stability 
of  the  model  as  the  cone  and  canopy  are  moved  to  varying  H/D^^ 
positions,  as  defined  in  Pig  7-5  of  Progress  Report  No  20. 

The  effect  of  cone  size  is  investigated  by  using  two 
different  cone  sizes  for  each  test  position.  Bach  cone  has  a 
half  angle  0  of  34°.  For  the  "small"  cone  the  base  diameter  d 
is  37#  of  the  canopy  Inlet  diameter  D^  and  the  "large"  cone  has 
a  base  diameter  that  Is  47#  of  the  inlet  diameter. 

4.4.3  Results 

The  criteria  for  the  degree  of  stability  in  the  rigid 
model  tests  is  a  combination  of  model  vibration  and  the  shock 
pattern  stability  as  determined  from  analysis  of  the  Schlieren 
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movies.  The  results  of  this  analysis  are  presented  in  Table  7-1  • 
As  different  model  configurations  were  tested  under  this  project, 
the  following  trends,  first  noted  in  Progress  Report  No  20,  held 
trues 

a)  The  overall  stability  of  the  flow  patterns 
increased  with  increasing  geometric  porosity  of 
the  models. 

b)  As  the  distance  from  the  cone  to  the  canopy  was 
increased,  stability  of  the  flow  pattern  decreased. 

For  a  third  trend  that  was  observed  for  models  of  lower  geometric 
porosity,  the  40#  model  displayed  a  contrary  effect.  In  this 
case,  the  model  was  more  stable  with  a  large  cone  than  with  a 
small  cone.  It  was  also  noted  that  the  “buzzing"  effect  repre¬ 
sentative  of  lower  porosity  models  was  absent  in  tests  on  the 
40#  model,  indicating  that  mass  was  not  being  spilled  over  the 
front  of  the  canopy  as  it  was  in  the  models  with  lower  geometric 
porosity. 

4.4,4  Proposed  Work 

In  view  of  the  above  experiments,  it  is  proposed  that 
rigid  pressure  distribution  models  of  the  40#  geometric  porosity 
configuration  be  constructed.  These  models  will  then  be 
tested  at  the  H/D^  positions  that  have  shown  the  best  results 
to  date.  Experiments  will  also  be  conducted  to  investigate 
the  characteristics  of  higher  geometric  porosity  rigid  models. 
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BUTTERFLY  STITCHING 


FIG  7-1  DETAIL  VIEW  OF  ATTACHMENT  OF  SUSPENSION  LINES  TO 
CANOPY  OF  SPIKED  PARACHUTE 


PIG  7-2.  FLEXIBLE  SPIKED  SOLID  CLOTH  CANOPY  PARACHUTE 
AT  MACH  3.0. 

{H/Di  =  0.62,  Do/Di  =  0.90,  l/D±  =  O.35) 
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FIG  7-3.  WIND  TUNNEL  SUSPENSION  ARRANGEMENT  FOR 
FLEXIBLE  SUPERSONIC  PARACHUTE  MODEL  DRAG 
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FIG  74.  WIND  TUNNEL  TEST  ARRANGEMENT  FOR 

FLEXIBLE  SUPERSONIC  PARACHUTE  MODEL 
WITH  FOREBODY 


FIO  7-5.  BALLISTICS  RANGE  TEST  DEVICE  FOR  RIGID  MODEL  OF  SPIKED  SOLID  CANOPY  PARACHUTE 
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5.0  Theoretical  Analysis  of  the  Dynamics  of  the 

Opening  Parachute 

No  work  concerning  this  project  will  be  reported 
during  the  present  reporting  period. 
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Project  No  9 


6.0  Statistical  Analysis  of  Extraction  Time,  Deployment 

Time,  Opening  Time,  and  Drag  Coefficient  for 
Aerial  Delivery  Parachutes  and  Systems 

A  draft  of  the  technical  report  entitled,  "Analysis 
of  Parachute  Deployment  of  Heavy  Cargo"  by  Helmut  G.  Heinrich, 
Ivan  R.  Russell,  and  Donald  J.  Eckstrom  has  been  submitted  to 
the  Procuring  Agency  for  review. 
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7.0  Gliding  Aerodynamic  Decelerator 

7.1  Introduction 

The  objective  of  this  study  Is  to  develop  a  self- 
inflating  aerodynamic  decelerator  with  a  lift  to  drag  ratio 
of  two.  With  such  a  lift  to  drag  ratio,  a  parachute  would 
glide  at  an  angle  of  26.5°  with  the  horizon  and  would  have  an 
angle  of  attack  OC  of  63.5°.  Figure  12-1  shows  the  forces  and 
coordinates  for  a  gilding  parachute.  Progress  Report  No  20 
presented  an  account  of  experimental  investigation  of  gliding 
characteristics  of  present  aerodynamlcal<iy  unstable  parachutes 
and  conventional  steerable  parachutes.  None  of  the  parachutes 
tested  met  the  necessary  glide  angle  specifications. 

Experimental  balsa  wood  and  clay  models  of  new  con¬ 
figurations  constructed  during  previous  reporting  periods 
achieved  a  glide  angle  close  to  that  desired.  These  models  as 
noted  in  Progress  Reports  Nos  18  and  19 »  were  improved  upon 
and  during  the  last  reporting  period  a  specific  configuration 
was  described.  Drawings  of  the  solid  test  models  are  shown  in 
Figs  12-8  and  12-9  of  Progress  Report  No  20. 

1*2  Continuation 

t 

A  wood  and  metal  form  based  on  a  modified  37.5  injch 
nominal  diameter,  10#  extended  skirt  parachute  canopy  was 
constructed  and  is  shown  in  Fig  12-2.  With  this  model  it  was 
possible  to  obtain  gore  and  rib  patterns  of  the  proposed  gliding 
parachute . 

During  the  present  reporting  period  the  gore  and  rib 
patterns  were  made  and  a  cloth  parachute  was  constructed  from 
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1.1  oz  nylon  cloth  with  a  nominal  porosity  of  120  ft^/(f t^-min) . 
The  design  Is  based  on  a  10#  extended  skirt  canopy  modified  with 
a  gliding  surface  as  shown  in  Pigs  12-3  and  12-4.  The  inflated 
parachute  in  the  wind  tunnel  is  shown  in  Pig  12-3. 

7.3  Tests 

Tests  of  the  cloth  gliding  parachute  were  conducted 
in  the  open  section  of  the  subsonic  wind  tunnel  (Fig  12-6). 

The  model  was  tested  at  a  flow  velocity  of  approximately  25 
ft/sec.  A  special  suspension  system  was  developed  which  allows 
for  the  systematic  inflight  adjustment  of  individual  suspension 
line  length. 

The  wind  tunnel  tests  showed  the  design  of  the  orig¬ 
inal  configuration  to  be  basically  sound.  By  adjustment  of  the 
suspension  line  lengths  a  stable  angle  of  attack  CCof  25°  +3° 
was  achieved,  corresponding  to  a  glide  angle  with  the  horizon 
of  65°  +  3°  and  a  lift  to  drag  ratio  of  0.47.  To  achieve  test 
results  closer  to  the  desired  glide  characteristics,  this 
original  test  model  was  modified. 

A  three  inch  wide  strip  of  cloth  was  removed  from 
the  base  of  eight  of  the  gores  contained  in  the  10#  extended 
skirt  portion  of  the  canopy.  Slots  were  cut  from  four  of  the 
guide  surface  panels  along  the  line  of  intersection  of  the  roof 
and  guide  surface  panels.  The  modifications  allow  for  greater 
inlet  area  at  the  base  of  the  extended  skirt  canopy  and  provide 
an  exhaust  jet  area  at  the  rear  of  the  glide  surface. 

The  same  suspension  line  system  and  test  arrangement 
were  used  in  wind  tunnel  tests  of  the  modified  parachute  model. 
Pig  12-7  shows  the  modified  model  in  the  wind  tunnel.  The 
maximum  stable  angle  of  attack  was  generally  greater  than  that 
of  the  original  model.  By  shortening  the  forward  suspension 
lines  by  0.18  where  B0  is  the  original  line  length,  a 
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maximum  angle  of  attack  of  35°  +  3°  was  obtained  giving  a  glide 
angle  with  the  horizon  of  35°  ±3°.  In  this  attitude  the  para¬ 
chute  was  moderately  stable.  At  angles  of  attack  greater  than 
35°  the  front  portion  of  the  canopy  tended  to  deform  and  at  an 
angle  of  attack  greater  than  45°  the  canopy  collapsed. 

Further  tests  were  conducted  on  the  modified  model  to 
determine  what  effect  the  porosity  of  the  canopy  h^d  on  test 
results.  In  order  that  the  same  modified  model  could  be  used 
in  these  tests,  the  nominal  porosity  was  reduced  by  starching 
the  canopy.  As  is  customary,  the  low  porosity  model  displayed 
less  stability  than  did  the  more  porous  model.  In  tests *  the 
less  porous  starched  model  had  a  glide  angle  with  the  horizon 
of  30°  +  3°. 

7.4  Proposed  Work 

Further  modification  and  tests  will  be  made  on  the 
nylon  cloth  parachute  during  the  next  reporting  period.  In 
addition,  an  investigation  will  be  conducted  to  determine  the 
relevant  characteristics  of  the  A/P  28  S-3  (Air  Force  Systems 
Command  publication  AFFTC-TR-61-37,  August,  196l)  "HALO"  para¬ 
chute  which  incorporates  a  modified  MC-1  canopy  with  a  single 
orifice.  Also,  other  existing  parachutes  will  be  studied  when 
models  and  other  pertinent  information  become  available. 
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FIG  12-1  FORCES  AND  COORDINATES  FOR  A 
GLIDING  PARACHUTE 


PIG  12-2.  MODIFIED  37.5  INCH  NOMINAL  DIAMETER  10# 
EXTENDED  SKIRT  CANOPY  FOR  OBTAINING 
GORE  AND  RIB  PATTERN  OF  GLIDING 
PARACHUTE 


GLIDING  SURFACE 


■1/2-10%  EXTENDEI 
SKIRT  CANOPY 


ROOF  PANELS 


GUIDE  SURFACE 
PANELS 


FIG  12-3. GLIDING  PARACHUTE  (SIDE  VIEW) 


RIB  -  6 


RIB-4 


RIB  -3 


RIB-2 


RIB  - 


FIG  12-4*  GLIDING  PARACHUTE  (TOP  VIE 


FIG  12-5  .  GLIDING  AERODYNAMIC  DECELERATOR  IN  WIN'D  TUNNEL 
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8,0  Effective  Porosity  Studies 

8.1  Parachute  Cloths 

Progress  Report  No  20,  Section  9.1»  presented  results 
obtained  from  tests  conducted  on  specially  mounted  samples  of 
four  different  parachute  cloths,  relating  effective  porosity,  C, 
to  the  pressure  ratio,  4  P/kPcr  *  with  density  ratio.  O',  as  a 
parameter.  The  special  mounting  was  used  to  eliminate  stretch¬ 
ing  and  slipping  of  the  cloth  sample  in  its  test  frame.  Despite 
the  special  mountings,  previously  established  trends  were  followed 
and  these  tests  served  merely  as  a  validation  of  earlier  experi¬ 
ments  . 

During  this  reporting  period  no  new  results  are  pre¬ 
sented.  Further  research  will  remain  on  a  relatively  low 
priority  basis. 

8.2  Ribbon  Grid  Configurations 

8.2.1  Introduction 

During  the  last  reporting  period,  tests  relating 
effective  porosity  to  pressure  ratio  with  the  grid  angle  0 
as  a  parameter,  were  conducted  on  a  series  of  ribbon  grid 
patterns  fabricated  from  steel  shim  stock.  Through  the  use 
of  this  relatively  rigid  material,  some  of  the  factors 
associated  with  cloth  deformation  were  eliminated.  The  tests 
were  conducted  on  grids  with  included  grid  angles  0  of  0°, 

15 °»  45 60 °,  75 °»  and  90°.  Aside  from  angle  variation,  all 
models  were  identical  in  construction,  each  having  a  geometric 
porosity  of  25#.  The  design  of  a  representative  grid  pattern 
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is  shown  in  Progress  Report  No  20,  Pig  13-6. 

During  the  present  reporting  period  the  data  from 
previously  conducted  effective  porosity  tests  of  1/4"  MIL-R- 
5608B,  Class  B  nylon  ribbon  grids  has  been  reduced.  By 
comparing  these  nylon  ribbon  grid  results  with  the  results 
from  steel  ribbon  grid  tests,  conclusions  can  be  drawn  concern¬ 
ing  the  effect  of  cloth  deformation  on  effective  porosity. 

8.2.2  Results 

The  data  from  the  tests  of  steel  and  nylon  ribbon 
grids  is  presented  in  Pigs  13-1  and  13-2.  The  effective  porosity 
as  a  function  of  the  pressure  ratio  with  the  grid  angle  ©  as  a 
parameter  varies  in  a  similar  manner  for  both  the  nylon  grids 
(Fig  13-2)  and  the  steel  grids  (Pig  13-l)  for  O.O34:  A  p/apcr^ 
2.00,  with  the  effective  porosity  being  generally  greater  for 
the  nylon  grids . 

In  Pigs  13-3  and  13-4  the  same  data  as  in  Pigs  13-I 
and  13-2  is  plotted,  but  with  the  effective  porosity  as  a 
function  of  the  grid  angle  ©,  Ap/apcr  being  the  parameter. 

In  these  two  graphs  the  behavior  of  the  curves  for  the  nylon 
and  for  the  steel  grids  shows  a  marked  difference.  In  this 
case,  for  0 0  £  Q  $  90°,  the  effective  porosity  of  the  steel 
ribbon  grids  increases  steadily  with  increasing  angle  0  (Fig 
13-3) •  Since  all  of  the  grids  have  a  geometric  porosity  of 
25#  this  variation  Indicates  that  the  discharge  coefficient 
of  the  grid  openings  increases  as  ©  Increases .  A  difference  is 
seen  in  the  behavior  cf  the  effective  porosity  versus  grid 
angle  curve  for  the  nylon  grids  (Fig  13-4).  For  45°  ^  ©  ^  90° 
the  effective  porosity  increases  as  ©  increases,  while  for 
0°  4  6  4  45°  the  effective  porosity  is  larger  for  smaller 
angles . 

This  latter  characteristic  is  explained  by 
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investigation  of  the  deformation  of  the  nylon  grids  during 
testing.  Figure  13-5  shows  photographically  this  deformation 
for  30°,  43°,  and  90°  angle  nylon  grids.  For  nylon  ribbon 
grids  with  0  4 30°  there  is  a  particularly  large  deformation 
of  the  grid.  This  results  in  a  net  increase  in  geometric 
porosity,  thereby  explaining  the  variation  shown  in  the  graphs. 
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FIG  13-1.  EFFECTIVE  POROSITY  OF  STEEL  RIBBON  GRID 
CONFIGURATIONS  AS  A  FUNCTION  OF  DIFFERENTIAL  PRESSURE 


FIG  13-2.  EFFECTIVE  POROSITY  OF  NYLON  RIBBON  GRID  CONFIGURATIONS 
AS  A  FUNCTION  OF  DIFFERENTIAL  PRESSURE 


©  =  30° 


PIG  13-5.  DEFORMATION  OF  NYLON  RIBBON  GRIDS 

AT  ap/a  Pcr  «  1.0 
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9.0  Study  of  Flow  Patterns  of  Aerodynamic  Decelerators 

»ty  Means  of  the  Surface  Wave  Analogy 

Izl  Introduction 

Progress  Report  No  20  reported  the  final  alignment  of 
the  water  channel  and  guide  rails,  the  completion  of  a  shadow¬ 
graph  and  a  Schlieren  system  for  detailed  flow  observation,  work 
concerning  the  design  of  special  water  depth  probes,  and  con¬ 
struction  of  a  model  of  the  supersonic  parachute  configuration 
proposed  in  Project  No  7.  Also  described  was  a  special  dye 
technique  for  obtaining  the  water  depth  profile  around  a  model. 

9.2  Present  Work 

9.2.1  Test  Section  Depth  Survey 

Before  conducting  any  quantitative  tests,  a  depth 
survey  of  the  channel  test  section  area  was  carried  out  to  deter¬ 
mine  the  tolerance  margins  of  depth  measurements.  For  this 
purpose  the  depth  probes  described  in  the  last  progress  report 
were  used.  With  the  channel  floor  set  horizontally  in  both  the 
longitudinal  and  transverse  directions,  the  water  depth  in  the 
test  section  area  was  found  to  be  within  +  0.003  inch  from  the 
reference  value. 

9.2.2  Channel  Slope 

The  flow  rate  in  the  channel  is  a  function  of  the 
water  depth  and  channel  slope  in  the  flow  direction.  The 
commonly  accepted  Manning's  formula  for  open  channel  flow 
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^Ref  l)  was  employed  to  determine  the  channel  slope  necessary 
to  produce  a  given  Proude  number  (Pr).  Hence  the  velocity 
corresponding  to  different  water  depths  is  determined  by. 


V  = 


ua  „„2/y/2 


(10.1) 


where  V  ~  water  velocity  in  feet  per  second 

Rh  =  hydraulic  radius  in  feet;  this  is  the  ratio  of 

flow  cross-sectional  area  to  the  wetted  perimeter 
S  =  slope  of  the  channel  bed 
n  =  roughness  coefficient. 

Prom  Ref  2,  a  value  of  n  =  0.008  was  used  for  the  smooth  glass. 

In  the  hydraulic  flow,  the  Froude  number  (Fr)  simulates 
the  Mach  number  (P)  of  the  aerodynamic  flow.  The  quantitative 
tests  reported  here  were  for  a  Proude  number  Fr  =  2. 


9.2.3  Instrumentation 

A  special  nozzle  depth  gage  to  enable  a  quick  setting 
of  the  opening  of  the  variable  nozzle  to  the  desired  free  stream 
water  depth  was  constructed.  It  consists  of  a  plexiglass  wedge 
with  a  rise  of  1  inch  over  a  length  of  10  inches  and  is  graduated 
along  the  length  in  0.1  inch  divisions , corresponding  to  a  0.01 
inch  rise. 

It  was  observed  that  the  plastic  Fresnel  lenses  of  the 
Shadowgraph  and  Schlleren  systems  tended  to  warp  In  their  orig¬ 
inal  mounting  rings  resulting  in  minor  distortions  and  a  poor 
field  image.  This  defect  was  remedied  by  remounting  each  Fresnel 
lens  between  two  sheets  of  single  weight  glass. 

A  large  piece  of  clear  plastic,  covering  the  field  of 
view  of  the  14  inch  diameter  lens, was  ruled  off  in  1/2  inch 
squares  and  attached  to  the  glass  under  the  test  section  area, 
thereby  providing  a  reference  grid  which  aids  in  the  Initial 
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positioning  of  the  model,  defines  clearly  the  probe  positions  in 
the  flow  field,  and  helps  interpret  the  shadowgraph  and  Schlieren 
images . 

Some  preliminary  tests  conducted  during  the  previous 
reporting  period  with  a  water  pressure  probe  consisting  essen¬ 
tially  of  an  impact  or  pitot  tube  connected  to  a  micromanometer 
showed  the  necessity  for  minimizing  the  size  of  the  probe  and  a 
careful  choice  of  its  dimensions.  A  recently  developed  pressure 
sensitive  paint  (Celab  Pressure  Sensitive  Paint  #  9-A)  appeared 
to  hold  great  promise  for  use  as  a  simple  transducing  element 
for  hydrostatic  and  total  pressure;  probes.  Preliminary  tests 
with  this  paint  were  carried  out  to  evaluate  its  general  operat¬ 
ing  characteristics  such  as  sensitivity,  repeatability,  and 
drift.  Qages  have  been  constructed  on  pieces  of  plastic  0.015 
inch  thick.  Unfortunately,  all  gages  constructed  so  far  have 
shown  a  marked  zero  di^ift  and  it  appears  that  the  design  of  a 
satisfactory  pressure  transducer  of  this  type  for  water  channel 
application  may  require  rather  extensive  development  work. 

9.2.4  Calculation  of  Simulated  Pressure  Distribution  from 

Water  Depth  Measurements 

9. 2.4.1  Summary  of  Analogical  Quantities 

The  theoretical  basis  of  the  analogy  and  the  assumptions 
underlying  it  are  well  known  and  have  been  presented  in  Refs  2 
and  3.  It  is  convenient  here  to  summarize  the  analytical 
quantities: 


Qas  Flow 

Two  dimensional, 

gas  with  £r  - 
Cy 

Velocity  ratio. 


compressible 

2 


Water  Flow 

Liquid  flow  with  free  surface 
in  gravity  field 


Velocity  ratio. 


v2 


w 
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Density  ratio. 


Water  depth  ratio. 


T2 

Temperature  ratio,  jjr- 


Pressure  ratio. 


Water  depth  ratio. 


£2 

dl 

d2 


(Water  depth  ratio) 


2 


2 


Velocity  of  sound,  a 


Velocity  of  surface  wave  propa¬ 
gation,  c  =Vsd 


Vo  v 

Mach  number,  M  =  —  Froude  number,  Fr  =  — 

a  c 

Shock  wave  Hydraulic  jump. 

For  flows  involving  shocks,  the  analogy  is  imperfect  for  two 
main  reasons: 

a)  For  air,  jf'  =  1.4,  whereas  the  analogy  assumes  a 
hypothetical  value  of  =  2 

b)  The  change  in  internal  energy  for  air  is  not 
equivalent  to  that  for  water. 

In  a  simple  quantitative  application  of  the  analogy, 
one  may  derive  the  air  density,  temperature  and  pressure  dis¬ 
tributions  directly  from  the  water  depth  distribution  using  in 
each  case  the  correspondence  Indicated.  It  will  then  be  found 
that  varying  degrees  of  agreement  will  result  and  that  the  best 
agreement  obtained  is  that  between  the  depth  ratio  and  the 
density  ratio. 

9. 2,4,2  Modification  of  the  Direct  Analogy 

A  modification  of  the  direct  analogy,  originally 
proposed  in  Ref  2,  produces  greatly  Improved  numerical  results. 
The  basis  of  this  modification  Is  the  use  of  the  depth  ratio 
to  determine  the  analogous  density  ratio  only.  All  the  other 
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necessary  aerodynamic  quantities  are  calculated  from  the  density 
ratio  by  means  of  well  known  aerodynamic  relations.  In  this  way, 
the  good  agreement  between  air  density  ratio  and  water  depth  ratio 
is  more  fully  exploited. 

In  all  the  numerical  applications  reported  here,  the 
modified  analogy  was  used. 

9. 2. 4. 3  Experimental  Procedure 

Since  no  direct  water  velocity  measurements  can  yet 
be  made, owing  to  the  lack  of  suitable  total  pressure  probes, 
the  simulated  flow  Mach  number  was  obtained  by  means  of  wave 
angle  measurements  using  a  symmetrical  diamond  airfoil  with  a 
wedge  angle  of  18°  at  both  leading  and  trailing  edges. 

The  channel  flow  was  adjusted  to  the  necessary  condi¬ 
tions  for  simulating  a  Mach  2  flow  by  adjusting  the  nozzle  flap 
to  the  desired  opening  height  by  means  of  the  nozzle  depth  gage, 
setting  the  longitudinal  channel  slope  to  the  required  value  as 
calculated  from  Eqn  10.1  ,  and  by  adjusting  the  pump  discharge 
by  means  of  the  Walworth  gate  valve  until  the  oblique  shock  angle 
coming  off  the  leading  edge  of  the  diamond  airfoil  corresponded 
to  the  value  for  Mach  2. 

The  model  was  then  placed  in  the  test  section,  carer 
fully  aligned  with  respect  to  the  reference  grid,  and  water 
depth  measurements  were  taken  at  all  -grid  points  in  the  area  of 
interest  and  also  at  points  of  intersection  of  the  reference  grid 
and  the  model  surface  as  close  to  the  model  as  the  probes  allowed. 
In  order  to  minimize  individual  differences  of  interpretation, 
all  depth  readings  of  a  given  test  were  made  by  the  same  opera¬ 
tor.  However,  tests  of  the  same  model  by  different  individuals 
produced  substantially  the  same  results.  The  measurements  of 
water  depth  on  the  model  itself  were  obtained  by  means  of  the 
Eosin  dye  technique  described  in  Progress  Report  No  20,  and. 
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in  this  case,  a  correction  factor  for  the  meniscus  was  deter¬ 
mined  experimentally  using  the  same  model  with  placid  water  of 
known  depth.  The  correction  factor  thus  established  reduced 
the  indicated  height  by  0.08  inch  for  the  flat  or  gently  curved 
surfaces  and  0.04  inch  ft>r  the  comers. 


9. 2. 4. 4  Data  Reduction  Procedure 


1)  The  values  of  Pii/P(p,  and  q* /P«j»  correspond¬ 

ing  to  the  desired  flow  Mach  number  are  taken  from 
standard  compressible  flow  tables  such  as  Ref  4. 
Subscripts co  and  T  indicate  the  free  stream  and 
total  conditions  respectively. 

2)  From  the  experimental  value  of  the  local  depth 

d  and  the  free  stream  depth  d»  ,  the  dimensionless 
water  deptn  ratio  d/d«*  is  calculated  and  inter¬ 
preted  as  the  air  density  ratio 

3)  The  value  of  is  calculated  by  means  of  the 

simple  relation.; 


e  /&.  -  do. 2) 


4)  Using  Ref  4,  any  aerodynamic  property  correspond¬ 
ing  to^/Pr  such  as  P/Prp,  T/Tj.  or  local  Mach  number 
can  be  obtained. 

5)  The  local  pressure  coefficient  can  then  be 
obtained  by  means  of  the  relation: 


(10.3) 


9.2. 5  Experimental  Results 
9.2. 5. 1  Diamond  Airfoil  Tests 
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The  main  object  of  testing  this  model  was  to  develop 
the  experimental  and  data  reduction  procedures  and  to  compare 
the  experimental  results  with  those  of  accepted  aerodynamic 
theories . 

Table  14-1  gives  the  experimental  data  and  calculated 
results  for  the  diamond  airfoil  model  at  Mach  2. 

Figure  14-1  provides  a  comparison  between  the  pressure 
coefficients  obtained  experimentally  and  the  values  calculated 
theoretically.  There  is  a  satisfactory  agreement  except  at  the 
corners  of  the  model  where  the  analogy  indicates  a  rounded,  more 
gradual  change  of  pressure  instead  of  the  abrupt,  step  function 
change  predicted  by  theory. 

9. 2. 5. 2  Test  with  Model  of  Proposed  Supersonic  Parachute 

The  main  geometric  characteristics  of  this  model  are 
given  in  Fig  14-2A  and  Illustrated  in  the  accompanying  photographic 
view  (Fig  14-2B). 

Depth  probings  were  taken  at  all  the  reference  grid 
points  and  also  at  points  close  to  the  model  surface  both  inside 
and  outside  the  canopy.  The  experimental  and  data  reduction 
procedures  outlined  previously  were  followed  and  the  pressure 
coefficients  were  calculated.  Shadowgraph  and  Schlieren  pictures 
were  taken.  Figure  14-3  shows  the  point  designation  system  and 
the  calculated  contour  lines  of  constant  pressure  (isobars). 

Figure  14-4  is  a  shadowgraph  of  the  flow  on  which  the  contour 
lines  of  Fig  14-3  have  been  superimposed.  Figure  14-5  is  a 
Schlieren  image  of  the  same  flow  with  the  calculated  constant 
pressure  lines  also  superimposed.  It  is  interesting  to  note  the 
generally  good  correspondence  between  the  contour  lines  and  the 
Schlieren  image.  It  should  be  noted  that  Fig  14-5  was  obtained 
by  enlargement  from  a  small  part  of  a  35  mm  negative;  a  relatively 
long  exposure  time  of  1/2  sec  was  used  to  minimize  high  fre¬ 
quency  random  changes.  With  modifications  in  the  photographic 
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technique  and  simple  optical  refinements,  it  is  anticipated 
that  much  improved  shadowgraph  and  Schlieren  images  will  be 
obtained. 

9.2.6  Other  Experiments  with  the  Analogy 

9. 2.6. 1  The  Starting  FTow  of  the  Supersonic  Parachute 

While  investigating  the  supersonic  parachute  configura¬ 
tion  by  means  of  the  Schlieren  system,  certain  interesting  flow 
characteristics  were  observed  during  the  starting  phase;  they 
are  illustrated  in  Figs  14-6  A,  B,  and  C.  When  the  water  flow 
initially  passes  over  the  symmetrical  wedge  forming  the  spike, 
a  typical  wave  pattern  consisting  of  an  oblique  shock  at  the 
nose  and  a  Prandtl -Meyer  expansion  fan  at  each  wedge  corner 
is  established.  The  presence  of  the  canopy  creates  a  back  pres¬ 
sure  behind  the  spike  and  generates  a  compression  wave  which 
travels  upstream.  This  wave  is  clearly  shown  in  Fig  14-6A.  As 
this  wave  travels  upstream,  it  acts  on  the  expansion  waves  com¬ 
ing  off  the  wedge  corners  (Fig  14-6B),  eventually  eliminating 
these  expansion  waves  and,  in  their  place,  oblique  shocks  develop; 
in  stable  configurations,  these  waves  pass  ahead  of  the  canopy 
leading  edges  (Fig  14-6  C).  These  oblique  shocks  may  be  con¬ 
sidered  as  resulting  from  the  stagnant  wake  behind  the  wedge 
forming  a  virtual  extension  of  the  wedge  profile. 

9. 2. 6. 2  Effect  of  a  Primary  Body  on  the  Performance  of  the 

Supersonic  Parachute 

It  is  well  known  that,  under  certain  operating 
conditions,  the  presence  of  a  primary  body  may  have  a  marked 
influence  on  the  performance  of  trailing  deceleration  devices. 
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During  the  present  reporting  period,  a  preliminary  investigation 
of  the  effect  of  primary  body  wake  on  the  supersonic  parachute 
was  initiated.  For  this  purpose,  a  two  dimensional  model  of  an 
ogive  cylinder  wa3  used.  The  main  geometric  characteristics  of 
the  configuration  tested  are  shown  in  Fig  14-7A.  Figure  14-7^ 
is  a  Schlieren  photograph  of  the  wave  pattern  of  the  supersonic 
parachute  model  at  M  =  2  without  a  primary  body.  Figure  14-7C 
shows  the  wave  pattern  of  the  same  model  behind  a  primary  body 
with  D^/Dg  =  11/16  and  L/D^  =5.3.  It  is  apparent  in  this  case 
that  the  presence  of  the  primary  body  wake  causes  detachment  of 
the  shock  wave  at  the  spike  nose  and  creates  irregularities  in 
the  flow  with  a  consequent  loss  of  stability. 

9.3  Proposed  Work 

The  usefulness  and  versatility  of  the  newly  developed 
stationary  model  water  channel  facility  is  adequately  illustrated 
by  the  examples  presented  in  this  reporting  period,  particularly 
when  considering  the  quantitative  pressure  distribution  experi¬ 
ments.  Further  development  and  refinement  of  instrumentation 
and  experimental  techniques  will  be  introduced  whenever  the 
need  arises.  A  new  light  source,  U3ing  a  carbon  arc  has  been 
acquired  and  is  expected  to  greatly  improve  the  shadowgraph 
and  Schlieren  images,  allowing  for  movie  camera  recording  of 
unsteady  phenomena.  Attention  will  be  devoted  to  the  design 
and  construction  of  miniature  total  pressure  probes  to  enable 
detailed  velocity  surveys  and  complete  channel  calibration. 

In  conjunction  with  such  developments  and  refinements, 
further  investigations  with  the  proposed  supersonic  parachute 
configuration,  both  with  and  without  primary  bodies  will  be 
conducted. 

The  water  channel  as  presently  developed  is  also 
available  for  other  investigations  that- may  be  required  in 
connection  with  the  total  contract  effort. 
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TEST  CONDITIONS 

FROUDE  NUMBER:  Fr  =  2.0 
FREE  STREAM  WATER  DEPTH  :  d«  =  .304  in 
FROM  NACA  TR1135  FOR  M«*2 
8/ft*.230  P,/^=.128  <3,^  =.35  8 


TABLE  14-1.  EXPERIMENTAL  DATA  AND  RESULTS  FOR 

DIAMOND  AIRFOIL  MODEL 
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EFFICIENTS  WITH  THEORETICAL  VALUES 
R  A  SYMMETRICAL  DIAMOND  AIRFOIL  AT 


Dj  =3.57 
Db  =  185 
L/Dj  =0.35 


A)  GENERAL  ARRANGEMENT  AND  DIMENSIONS 


B)  PHOTOGRAPH  OF  WATER  CHANNEL  MODEL 


PIG  1 t-2.  SPIKED  SUPERSONIC  1 AR A  CHUTE  CONFIGURATION 


CONFIGURATION  AT 


A)  INITIAL  STAGE  -  NOTE  PRANDTL -MEYER  EXPANSION 
AT  CORNERS  OF  SF IKE  WEDGE 


D)  INTERMEDIATE  STAGE  -  EX  F  AMS  I  ON  FAN  STILL  FRESENT 


C)  FINAL  STAGE  -  STAPLE  FLOW  PATTERN  ESTABLISHED, 
EXPANSION  FAN  NOW  REPLACED  BY  SHOCK  WAVE 


FIG  !•.-  .  STARTIN',  FL  W  FOR  SPIKED  IAEA  CHUTE  CONFIGURATION  AT  M  -  2 


FLOW 


A)  GENERAL  ARRANGEMENT 


B)  FLOW  PATTERN  WITH  NO  PRIMARY  BODY 


C V  FLOW  PATTERN  WITH  PRIMARY  BODY:D1/D^,=  li/l6;L0/D«  5.3 


FIG  14-7.  EFFECT  OF  PRIMARY  BODY  ON  WAVE  PATTERN'  OF  SPIKED  PARACHUTE 


Project  No  16 


10.0  Stress  Analysis  of  the  T-10  Troop  Parachute 

Progress  Report  No  18  presented  an  analytical  method 
for  obtaining  the  stress  distribution  over  a  fully  inflated 
canopy  in  steady  descent. 

Efforts  have  now  been  directed  at  extending  the  basic 
theory  to  the  more  complex  case  of  finding  the  stress  over  an 
opening  canopy.  It  is  hoped  that  when  this  method  is  completely 
developed,  the  fully  opened  canopy  can  be  treated  as  a  special 
case  under  this  more  general  theory.  Certain  simplifications 
appear  possible  in  the  theory  as  thus  far  developed  and  these 
will  be  incorporated  in  the  general  development.  Preliminary 
results  of  these  new  efforts  should  be  available  during  the 
next  reporting  period. 
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